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ABSTRACT 

We analyse the one-dimensional anharmonic oscillator using effective operator meth- 
ods in both the strong and weak coupling limits. We show that in the case of a one- 
dimensional model space, the similarity transformation needed to define the effective 
Hamiltonian is related to the coefficients in the expansion of the wave function in the 
unperturbed harmonic oscillator basis. We obtain an infinite system of equations which 
is equivalent to those obtained from the Hill Determinant solution of the anharmonic 
oscillator. The analytic properties of the resulting equations reveal the non-perturbative 
features of the underlying problem. Thus, we demonstrate the the utility of the effective 
operator method for solving a non-analytic strong coupling problem. 
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The effective operator method [TJ has been used extensively and successfully within 
a cluster approximation scheme to obtain the low-lying spectroscopy of complex nuclei 
[2] with realistic nucleon-nucleon interactions. Central to the methodology of £Q is the 
iterative construction of a similarity transformation which transforms the original hamil- 
tonian to a new hamiltonian having a two-component block diagonal structure where 
one component is finite dimensional and accommodates the low-lying spectroscopy. Di- 
agonalising this finite dimensional block diagonal sub-matrix yields a finite number of 
eigenvalues to any desired precision corresponding to a subset of the exact solutions. In 
a sense, the similarity transformation is designed to decouple a finite dimensional sub- 
space from the rest of the spectrum, even in cases when the original hamiltonian cannot 
be treated by perturbative methods. However, there are few rigorous results on the ex- 
istence, or non-existence [3], of the similarity transformation on which the utility of the 
method hinges, although several practical issues specific to strongly correlated nuclear 
systems are under investigation. Thus the formal properties of this effective operator 
approach require additional study, particularly within the framework of problems known 
to be non-perturbative in character. In addition to formal aspects, a deeper understand- 
ing of the errors associated with various approximations is required. For example, when 
a sequence of clusters is introduced as in the nuclear many body applications ,2 a we still 
need to understand how to optimize the convergence with increasing cluster size and/or 
with increasing model space size. In order to clarify the physical utility of the approach 
and refine our knowledge of its limitations and properties we address a well studied 
non-perturbative problem, the one dimensional quartic oscillator in both the strong and 
weak coupling regime. In doing so, we shed light on the origins of the decoupling which 
plays such a crucial role. 

The hamiltonian H that we select, has the form H + V where 

2 2 2 

P . mu z x z 
2m 2 

and V = Xx 4 . 

The perturbation expansion for energy eigenvalues is known not to converge, inde- 
pendent of the size of A j3]. The divergent behaviour of the expansion may be traced 
to the fact that, to large orders in perturbation theory, the growth in the number of 
contributions is sufficiently rapid to lead to a series that eventually diverges, even for A 
arbitrarily small. 

Various non-perturbative methods have been applied to extract the eigenvalues [5]. 
In this investigation we will construct an effective Hamiltonian using a similarity trans- 
formation along the lines proposed by Lee & Suzuki pQ . In addition to extracting energy 
eigenvalues from the effective hamiltonian, the matrix elements of the generator of the 
similarity transformation will also be used to develop a system of equations equivalent 
to those obtained from the Hill determinant solution of the anharmonic oscillator. 

We begin by establishing the following notation: let E denote a generic eigenvalue 
and | \I/ > the corresponding exact eigenstate of the Hamiltonian in Eq.l. In that case 
we have 

H | * >= E | * > (2) 

Following |T], we define S, the generator of a similarity transformation such that H = 
e -s e s anc j \y > = e ~s I vl/ > . It then follows that 

H\ > = E\ * > (3) 
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For now, S is arbitrary but subsequently, restrictions will be placed on S. 

The next step in the treatment along the lines of pQ is the identification of a suitable 
model space (P-Space). We will use a one dimensional model space containing the 
ground state of H signified by | >. With this choice, P, the projection operator onto 
the model space is just | >< |. The Q space is then E | i >< i |, where the sum over 
% runs over all states of the harmonic oscillator other than the ground state; P and Q 
are orthogonal and P + Q = I . 

The first constraint we impose on S is that S = QSP. This leads to the identities 

PS = SQ = 

and 

S = SP = QS 

These identities will be used extensively later on. Since P and Q are orthogonal, S 2 = 
leading to e ±s — I ± S. With this choice of P and Q spaces and the restriction on 
S, S reduces to a column vector whose non-vanishing matrix elements are of the form 
< i | S | > (i / 0) which will be denoted by S^. The Si will be seen to be related to 
the expansion coefficients of the true wavefunction in the unperturbed basis. The final 
condition we impose on S comes from requiring that 

HP\ ^ > = EP\ ^ > (4) 

As in we define the effective hamiltonian H e ff to be PHP. In general, diagonal- 
ising H e ff yields a subset of energy eigenvalues of the full hamiltonian, since in our case 
H e ff is one dimensional we will recover just one eigenvalue for any given 5*. Expanding 
H e fj in terms of S gives 

H eff = PHP + PVS (5) 

From Eq. 0]and the fact that P and Q are projection operators it follows that QHP = 0, 
yielding 

SHP + SVS = QVP + QHS (6) 

which can be rewritten as 

QVP + QHS = SH eff 

After subtracting Q S from both sides of the preceeding equation, (Q is an arbitrary 
parameter), we obtain the following equation for S 

S = A(QVP-S(H eff -n)) (7) 

where 

A = (O - QHQY 1 

Note that since S and H e ff must be Q independent, Q is strictly speaking redundant. 
Nonetheless, the utility of Q will soon be made clear. No use has been made of the 
fact that the P space is one dimensional, thus Eqs. EJ El & El ar e perfectly general and 
identical equations arise when treating more complex hamiltonians along the lines of ^Q. 

Since H e ff depends explicitly on S, Eq. [7| is a non-linear operator equation for S 
potentially admitting more than one solution. In addition, A is ill-defined in an infinite 
dimensional Hilbert Space. In applications to many body problems it is customary 
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to impose an arbitrary cut-off in the Q space which leads to a well defined A, in our 
case this amounts to neglecting all oscillator quanta above a certain energy. This also 
reduces S to a finite dimensional column vector. Even with this simplification Eq. [7| is 
still hard to solve exactly; in applications to many body hamiltonians iterative solutions 
are preferred. We propose the following sequence of iterations; 

= A{QVP-S {n - l \H^ n \ ff -^l)) (8) 
H {n+l) eff = PHP + PVS {n) (9) 

with = as the starting point in the sequence. After a suitable number of iterations 
the sequence is terminated. Since the P space is one dimensional, the final H e ff is simply 
the energy eigenvalue, no diagonalisation is necessary. On the other hand, it is necessary 
to check the cut-off independence of the energy eigenvalues and the convergence of the 
iterations, this can be done by varying both the cut-off and Q. As a check we calculate 
the eigenvalues emerging from this procedure after a total of seven iterations. We retain 
only the first ten harmonic oscillator eigenstates of even parity in calculating A; since 
V does not mix eigenstates of even and odd parity it is adequate to consider only even 
parity states in calculating A. As in jU] we use h = 1, u = 2 and m = .5 and allow A to 
vary. Our results (setting Q, — — 1) are 



Iteration(n) 


A = .1 


A = .2 


A = .3 


A = 1 


A = 2 


A = 3 





1.075 


1.15 


1.225 


1.75 


2.5 


3.25 


1 


1.06792 


1.12612 


1.17806 


1.45479 


1.7359 


1.95799 


2 


1.06603 


1.12035 


1.16752 


1.40423 


1.62662 


1.79234 


3 


1.06551 


1.11886 


1.16495 


1.39484 


1.61088 


1.77311 


4 


1.06535 


1.11845 


1.16429 


1.3929 


1.60819 


1.77022 


5 


1.06531 


1.11834 


1.16411 


1.39248 


1.60768 


1.76971 


6 


1.06529 


1.11831 


1.16407 


1.39238 


1.60757 


1.76962 



These results are in excellent agreement with the ground state energies in jB] ; extending 
the size of the Q space used to obtain A changes the eigenvalues (if at all) only at the 
fifth significant figure or beyond. Furthermore, allowing Q to vary between -2 to +2 
leads to a similarly small change in the obtained eigenvalues. As in [7] we can access 
different solutions for S and thus different eigenvalues, by choosing a different range of 
values for Q. 

Given that the hamiltonian is non-perturbative the convergence of the sequence of 
iterations is non-perturbative in character; furthermore the error induced by truncating 
the Q space to construct A is small. Since truncation of the Q space in the manner 
just described is a key ingredient in the application of effective operator methods in 
many body hamiltonians, we will attempt to understand why this procedure is valid. 
More precisely, we will show that at least for the purposes of calculating low lying 
eigenvalues, S may be accurately approximated by a column vector of finite dimension; 
since 5* = QSP, it then follows the Q space needed to define A is effectively finite 
dimensional which justifies the truncation we employ. 

We start by expanding Eq. 0] in terms of 5" and H leading to 

(H — SH + HS — SHS)P | >= EP | > (10) 
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In order to derive an identity for E in terms of the Si we expand | \& > in terms of 
the eigenstates of H as follows 

|^>=5>„|n> (11) 

The sum in Eq. begins from zero and runs over all positive integers. For the sake of 
future notational convenience we will denote < n | A \ m > by A nm where A is any 
operator. Inserting the form of | \? > from Eq. ^2 into Eq. and operating with < | 
on both sides yields 

#=<0|if|0> + <0|#S|0> (12) 
assuming that «o is non-vanishing. Inserting a complete set of states yields 

E = H 00 + V 02 S 2 + V04S4 (13) 

exploiting the fact that the quartic perturbation connects the ground state only to 
states I 2 > and | 4 >. Eq. EI] is a special case of a more general expression for the exact 
eigenenergy resulting from an effective hamiltonian in a one-dimensional model space 
derived in jjj. 

Equations for the remaining Si may be obtained by sandwiching Eq|U] between | > 
and < n | where n / 0. Doing so and using the fact that P and Q are projection 
operators gives 

<n \ SH I > + < n I SVS \ > = < n \ V \ > + < n \ HS \ 0> (14) 

By judicious insertions of complete sets of states, the left hand side of Eq. HUcan be 
reduced to 

S n {H 00 + V02S2 + V04S4) 
The term within brackets is just E from Eq. El Thus Eq. IT4lmay be written as 

S n E =<n|1/|0> + <n|iLS|0> (15) 

Inserting a complete set of states yields 

S n E = V n o + H nm S m (16) 

m 

As before we will assume n to be even, the extension to include odd n is straightforward. 
Inserting different values of n into Eq. HH1 leads to a system of coupled equations for the 
S n . For n = 2 we have, 

V 20 + (H 22 - E)S 2 + V 24 S 4 + V 26 S 6 = (17) 

For n = 4 we have, 

V A0 + V 42 S 2 + (#44 - £)S 4 + V A6 S 6 + V A8 S 8 = (18) 
For n = i and i > 4, there is no term independent of S. The equation takes the form 
Vii-4>S'j-4 + Vu-2Si- 2 + (Ha — E)Si + Vu +2 Si +2 + Va+^Si+i = (19) 
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We can make all the preceeding equations homogenous with the generic substitution 

c _ A 

With this substitution Equations ITBI IT71 [TH1 read 

(H o-E)p + V 02 p 2 + V 0i p i = (20) 
V 20 P + (H 22 -E)(3 2 + V 24 (3 i + V 26 (3 G = (21) 
V m p + V i2 p 2 + (H u -E)p 4 + V ie p 6 + V, 8 Ps = (22) 

For larger values of n it is adequate to replace Si by in Eq. HH] giving 

Vu-tPi-i + Vii_ 2 f3i_ 2 + (Ha — E)f3i + V ii+2 f3i +2 + V^ +4 /? i+4 = (23) 

For the sake of comparison, we substitute the expansion from Eq. ^2 into Eq. Re- 
stricting ourselves to a« with % even yields 

(H Q0 - E)a + V 02 a 2 + V 04 a 4 = (24) 
V 20 a + (H 22 - E)a 2 + V 2A a 4: + V 26 a e = (25) 
V m a + V A2 a 2 + (H u - .E)a: 4 + V m a & + V^ag = (26) 

and for larger values of i we have 

V^j_ 4 Q! i _4 + V^j_ 2 aj_2 + - E)ai + Vu +2 a i+2 + Vn +/L a i+ 4 = (27) 

It is clear that ai and satisfy the same set of equations. However, solving for the o.i is 
equivalent to solving the Hill determinant for the quartic oscillator j^j. Since the Si are 
proportional to the fy, the Si may also be obtained from solving the Hill determinant. 
Furthermore, the physical significance of the Si is apparent, the Si may be taken equal 
(up to an overall constant) to the coefficients arising in the expansion of the true wave 
function in the unperturbed basis. As far as we are aware, this is the first instance where 
a connection between S and the expansion of the wave function in unperturbed basis 
has been established. 

The ofj may be independently determined by the requirement that 

det(H — EI) = 

for an invertible H. It can be shown that this requirement is equivalent to the more 
general statement in [7j of which Eq.^Jis a special case [8 . This provides an independent 
cross-check of our results. 

There are additional constraints on the Si for large i which arise from the behaviour 
of the matrix elements at large i. To derive these constraints, let us rewrite Eq. EH 
as follows, 

Va-iSi-i + Vu^ 2 Si_ 2 + HaSi + Va +2 Si +2 + Vu^S^ = ESi (28) 
For large i the matrix elements appearing in Eq. EH have the approximate form 

< i | x 4 | + 4) > = (A)3(I)i2(i + 5 /i) + ... (29) 
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< % 


x 4 | 


(i + 


2) 


> = 


( ) 2 * 2 (1 + 3A) + --- 
moo 


(30) 




< % 


| x 4 


z 


> = 


( h ) 2 ( 3 )(2t 2 + 2t + l) 


(31) 


< i 


A i 

x 4 | 


(i- 


2) 


> = 


( * )V(1 + - 


(32) 


< i 


x A | 


(<- 


4) 


> = 


/loin 

— 2 t i 2 1-5/i +■■■ 
mcu 4 


(33) 


< i 


x 2 | 


(< + 


2) 


> = 


(34) 


< i 


x 2 \ 


(<- 


2) 


> = 


(^-)(2< + l) 


(35) 




< i 


|x 2 


i 


> = 


(36) 



where higher order finite terms in the expansion have been neglected. As can be seen 
from the above equations, the left hand side of Eq. |2H contains both quadratic and linear 
divergences in i at large i, and no such divergences appear on the right. Eq. |2Hl w ill be 
consistent only if Si fall off rapidly for large i aided by possible cancellations due to 
sign differences between the various Si appearing on the left. As a consequence of the 
Si falling off rapidly, low energy eigenvalues of the full hamiltonian are expected to be 
only weakly dependent on high energy eigenstates of the unperturbed hamiltonian. As 
a check, we allow i in Eq. |2H to run over a limited range of even values beginning at 
0, and compute the lowest energy eigenvalues from the resulting equations. We use the 
same numerical values for the parameters as before and allow A to vary over the same 
range. The results are summarised in the table below: 



imax/2 j 


A -> 


.1 


.2 


.3 


1 


2 


3 


5 




1.06529 


1.11829 


1.16406 


1.39337 


1.61123 


1.77481 


10 




1.06529 


1.11829 


1.16405 


1.39235 


1.60755 


1.76963 


15 




1.06529 


1.11829 


1.16405 


1.39235 


1.60754 


1.76959 


20 




1.06529 


1.11829 


1.16405 


1.39235 


1.60754 


1.76959 



The entries in the table are the lowest energy eigenvalues obtained from sets of not 
more than 20 equations, and are in excellent agreement with the results obtained by 
our iterative procedure and with [§]. This is strong support of our earlier claim that 
low energy eigenvalues of the full hamiltonian are only weakly dependent on high energy 
states of the unperturbed hamiltonian. It is now clear why the truncation of the Q space 
which was implemented in order to facilitate the iterative procedure is justified. The 
decoupling of high energy states in a toy model has been independently studied in [5]; 
however unlike the quartic oscillator the hamiltonian in jH] can be analytically inverted, 
permitting the use of methods very different from the ones we employ in this paper. 

As an additional analytical check on the decoupling of Si for large i, we consider the 
following Ansatz for large i, 

% ~ (-1)^(1 - MMrnW (37) 
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where / is some arbitrary unspecified function. Inserting this Ansatz into Eq. |2H1 and 
using the asymptotic forms of the matrix elements given earlier, we see that Eq. EHl is 
free of divergences in % provided 

/(-2) - /(2) + 2/(1) - 2/(-l) + 1 = (38) 

The utility of Eqs. EE & EH w iU become apparent when we investigate what happens 
when A is so large that it dominates the quadratic term in the potential. In that case 
the hamiltonian may be conveniently expressed as 

2 2 2 2 2 

p muj x 4 muo x 



Following the same procedure as before but with V = Ax 4 — mu) ^ x gives an equation 
identical to Eq. EHlbut with different values of Ha, and Ki±2- This is not surprising as 
the structure of Eq. |2H does not rely on the precise form of V but on the fact that the 
only non-vanishing matrix elements V^- have | (i — j) |< 4, which is the case for both 
forms of V that we consider. 

Once again, consistency requires that all divergences quadratic and linear in i cancel 
on the left hand side of Eq. EBJwith V modified to study the strong coupling limit. Using 
the approximate forms of the relevant matrix elements given earlier, it is straightforward 
to verify that the recursion relation in Eq. EH originally derived to analyse finite coupling 
is sufficient to gaurentee the cancellations in the strong coupling limit as well. The only 
difference between the two cases lies presumably only in the 0(l/i 2 ) term which plays no 
role in cancellation of divergences. This is a further indication of the non-perturbative 
nature of the effective operator method. 

To conclude, we have implemented an an effective operator treatment of the anhar- 
monic oscillator on the lines of PQ. We have analyzed the validity of a key simplification 
made in the treatment of the nuclear many body problem. The role of the expansion 
of the exact eigenfunction in the basis of the unperturbed hamiltonian in defining S, 
the generator of the similarity transformation has been emphasized. Both the strong 
and weak coupling cases may be treated along the same lines, underscoring the non- 
perturbative nature of the formulation. Our numerical application demonstrates rapid 
convergence to known results with modest effort. 

This work was supported in part by U.S. DOE Grant No. DE-FG-02-87ER-40371, 
Division of High Energy and Nuclear Physics. 
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